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Retinal implants for restoring partial vision to people 
with degenerative photoreceptor diseases consist of an 
array of electrodes implanted on the retina which can 
be driven according to images captured by a camera via 
stimulation electronics. In the past two decades this 
technology has progressed from prototypes, through 
clinical trials, through to the first commercial devices 
available in the U.S. last year (2013). Currently, the 
functional benefit to patients is largely as an aid to 
mobility and navigation. Improved outcomes, for exam- 
ple to the ability to perform useful object recognition, 
will likely require higher density electrode arrays. There 
is also likely to be benefit from attaining greater spatial 
control of electrically evoked neural activity by choosing 
appropriate stimulation parameters. As a first step 
towards this goal we investigated the ability of simple 
linear-nonlinear models to capture the response of 
visual cortical neurons driven by a retinal implant. 

Using an electrode array containing 21 platinum elec- 
trodes in the suprachoroidal space of the retina in an 
anaesthetized cat, spatially white and temporally sparse 
electrical pulses were delivered to the retina simulta- 
neously across several electrodes. Each pulse was either 
an anodic-first or cathodic-first charge-balanced biphasic 
current pulse. Multi-unit spiking responses were 
recorded using 6x10 arrays (Blackrock Microsystems, 
UT) in the primary visual cortex. 

We formulated a model to predict spiking responses 
observed in the cortex in response to electrical 
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stimulation of the retina. The pathway between the 
retina and the cortex was modelled as a set of spatial 
linear filters, each followed by a static nonlinearity. To 
encapsulate the effects of both anodic-first and catho- 
dic-first current pulses, two different linear filters (one 
for anodic-first and one for cathodic-first stimuli) and 
their associated nonlinearities acted on the input sepa- 
rately. The output of these two nonlinearities was com- 
bined to produce an estimate of the spike rate at a given 
site in the cortex. The parameters of the model were 
optimized with a non-linear fitting algorithm. We used 
83% of the available experimental data to fit the model 
parameters. After fitting the model parameters, the 
model accurately predicted the pulse-by-pulse spike 
count for the remaining 17% of the data (r=0.83). 

The spatial linear filter can be considered to be the 
weighting matrix for the effect of the stimulating elec- 
trodes on a cortical site. We anticipate that the model 
may provide a basis for devising ideal stimulation strate- 
gies for shaping the activity in the cortex. 
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